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We present a new method for extracting excited states from a single two-point correlation function 
calculated on the lattice. Our method simply combines the correlation function evaluated at different 
time slices so as to "subtract" the leading exponential decay (ground state) and to give access to 
the first excited state. The method is applied to a quenched lattice study (volume = 24^ x 64, 
(3 = 6.2, a~^ — 2.55 GeV) of the first excited state of the nucleon using the local interpolating 
operator O — eabclu^Cjsdbluc. The results are consistent with the identification of our extracted 
excited state with the Roper resonance A^'(1440). The switching of the level ordering with respect 
to the negative-parity partner of the nucleon, A'^*(1535), is not seen at the simulated quark masses 
and, basing on crude extrapolations, is tentatively expected to occur close to the physical point. 



1. INTRODUCTION 

Hadron spectroscopy offers a well defined framework 
to test the non-perturbative regime of QCD and to shed 
Ught on the mechanism responsible for confinement. Of 
course, the masses of ground states do not offer sufficient 
information to unravel the dynamics responsible for that 
mechanism, but fundamental insight can be gained by 
accessing the masses of excited states as well. 

As is well known, one of the long-standing puzzles in 
baryon spectroscopy is the level ordering between the 
negative-parity partner of the nucleon, iV*(1535), and 
the first excitation of the nucleon, the so-called Roper 
resonance iV'(1440). Phenomenological low-energy mod- 
els of QCD are generally quite successful in explaining 
many features of baryon mass spectra Q . Some of them 
are however unable to explain the proper level ordering 
between positive- and negative-parity excited states of 
the nucleon, like e.g. the one-gluon-exchange quark model 
of Ref . 2] ; such a failure has led to the speculation that 
the Roper resonance cannot be described by a picture 
based only on three constituent quarks. Other models do 
reproduce the correct level ordering, like the Goldstone- 
boson-exchange quark model of Ref. , where the level 
switching between the radial and orbital excitations of 
the nucleon is a consequence of the spontaneous break- 
ing of chiral symmetry. 

It is clear that a non-perturbative approach based on 
first principles, such as lattice QCD, offers a unique pos- 
sibility to unfold the above issues and, in particular, to 
investigate the evolution of the excitations of the nucleon 
as a function of the quark masses. Ground-state spec- 
troscopy in lattice QCD is well established and there is a 
general agreement with experimental data at the level of 
5-^ 10% precision 0]. On the contrary it is still challeng- 



ing to measure the mass of an excited state. Both ground 
and excited masses are extracted from the (Euclidean) 
time behavior of two-point correlators and therefore a 
clean separation among them is required. This problem 
is currently addressed following two main approaches^. 

The first approach uses a single two-point correlator, 
constructed from an interpolating field with appropriate 
quantum numbers. A naive multiple-exponential fit of 
the time dependence of the latter is in general affected 
by large systematic or statistical errors, due to the con- 
tamination from high excited states and/or to the large 
number of fitting parameters. To circumvent this prob- 
lem, methods inspired by Bayesian statistics, like the con- 
strained curve fitting method or the maximum entropy 
method 0, were developed and applied to lattice anal- 
yses aimed at studying the Roper resonance 8, 9]. The 
crucial point of these methods is clearly the choice of the 
priors (see Ref. lOj). 

The second approach uses several two-point correla- 
tors constructed from various interpolating operators Oi 
which couple differently to the states of interest Tl'l . A 
matrix of correlation functions is constructed choosing in 
all possible (independent) ways the source and the sink 
operators. By solving an appropriate eigenvalue problem, 
the contributions coming from the various states coupled 
to the operators Oi are disentangled, so that the masses 
of ground and excited states can be determined from the 
eigenvalues (see Refs. 0,^3)- The interpolating oper- 
ators must be chosen properly, taking into account the 
computational cost of constructing a matrix of correla- 
tors. The following two (local) operators are commonly 
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^ We mention that a throughout discussion of various mathemat- 
ical techniques for analyzing exponentially-damped time series 
can be found in Ref. 1^. 
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iployed miTlITi: 

O'(x) = 



Sabc[UaCj5db]Uc{x) 
= eabc[UaCdb]'y5Ucix) 



where C is the charge conjugation Dirac matrix. It is 
known that the first and the second operators in 
Eq. Hl.l|l are, respectively, strongly and very weakly cou- 
pled to the ground state of the nucleon. A recent im- 
provement [13 is the use of different smeared interpolat- 
ing operators as a tool to generate nodes in the radial 
wave function, which allows to capture more efficiently 
the Roper signal. The crucial point of the "matrix" ap- 
proach is a careful choice of the time interval, where the 
correlator matrix of dimension N has to be dominated 
by the first N states. 

We point out that a precise determination of excited 
states on the lattice is computationally very demanding 
on its own, since the size and the mass of the state in- 
creases with the excitation, which in turn implies the use 
of large lattice volumes and fine lattice spacings in order 
to keep under control systematic errors (see Ref. j^). 

In this paper we show that it is possible to access the 
first excited state coupled to a single two-point correlator 
G{t), defining a new correlator which is simply an appro- 
priate combination of G{t) taken at different time slices. 
The mass of the first excited state can be extracted from 
the plateau of the effective mass of the new correlator. 
Therefore our method allows to establish unambiguously 
the time interval where the ground and the first excited 
states dominate the two-point correlator G{t). In this 
way no prior is needed at all and the extracted resonance 
mass is not affected by contaminations from higher ex- 
cited states. 

Our procedure is explained in Section [51 while its ap- 
plication to the extraction of the Roper mass, using the 
local interpolating field O, is presented in Section O In 
Section^ we draw our conclusions and outlook. 



2. THE METHOD 

We limit ourselves to consider a single local opera- 
tor which has a good overlap with both positive- and 
negative-parity baryon states, i.e. 



0^{x) = SabclUaCj^db] Ucj{x) , 



(2.1) 



where Latin indices denote color and Greek indices re- 
fer to the Dirac structure, that is contracted in the (ud) 
diquark on the r.h.s. of Eq. (|2.1|) . Let us then define in 
the usual way the two-point (zero-momentum) correla- 
tion function for the operator H2.1|l as 

Gy'fit) = $^(0|T{Oy(x,i)a^(0,0)} |0) . (2.2) 



Adopting periodic boundary conditions on the lattice, 
Eq. 12.211 can be expanded as 



(1-1) G{t) = 




ZTe 



-aM+{T-t) 



-aM, t 



(2.3) 

where t is an integer number in units of the lattice spac- 
ing a, T represents the time extension of the lattice, the 
suffixes ± refer to the positive- and negative-parity states 
coupled to H2.1|l . respectively, Mf" refers to the masses 
of the given-parity tower of states, and Zf^ are the cor- 
responding couplings with the interpolating operator. 
Choosing 74 diagonal, one gets 

G{t) = ^[Gn{t) + G22{t) + G33{T - t) + G44{T - t)] 

(2.4) 



J2 [z+e-*^^ 



zre 



i>0 



Assuming (T — t) sufficiently large one has 

G(i)=5]z+e-*'^^ (2.5) 



from which the mass of the positive-parity ground state 
can be extracted from the large-time behavior through 
an effective mass analysis. 

Let us define the following quantity 



G{t) = G{t+l)G{t-l)-G^it) 



(2.6) 



which, by direct substitution of Eq. H2.5() . can be ex- 
panded as 



G{t) = 2 J2 Z+Z 



+ „-(aM++aM+)t 



[cosh(aM+ - aM+) - l] 



(2.7) 



One immediately sees that the diagonal terms i = j are 
cancelled out in the new correlator and that the lead- 
ing exponential decay in Eq. (|2.7|l is governed by the 
sum of the masses of the first two states qMq and aM^ . 
Thus, an effective mass analysis applied to G{t) provides 
the quantity a{M^ + M^). The use of the "modified" 
correlator 1)2. 7|l is expected to work when the mass gap 
between the ground state and the first excited state is 
sufficiently large, since for a small mass gap the factor 
[cosh(aMj^ — aMg) — 1], which multiplies the exponen- 
tial in the r.h.s. of Eq. H2.7|l . suppresses the signal. 

One can thus set up the following general procedure 
for extracting the first excited state from correlators 
like 
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i) Extract the mass of the ground state aM^ from the 
plateau of the effective mass for the correlator (|2.5|1 . 
Let the plateau region be given by the time interval 
{U,tf]. 

ii) Extract the sum of the masses of the first two 
states, a{MQ + M^), from the plateau of the ef- 
fective mass for the "modified" correlator (|2.6|l . Let 
the plateau region be the time interval [t* , i^] . This 
second plateau is always such that [t*, i^] < [ti, tf]. 
Indeed, in the time interval the correla- 
tor H2.5|l is dominated by the ground state sig- 
nal and the modified correlator (|2.6fl is pure noise, 
whereas in the time interval [t* , t^] the signal as- 
sociated to the first excited state is visible in both 
correlators. This in turn means that in the larger 
time interval [t* , t /] the correlator p. 5(1 is domi- 
nated by the first two terms due to the ground and 
the first excited states, while the contribution from 
higher excited states is negligible, i.e. within statis- 
tical fluctuations. 

iii) Check the consistency of the results obtained from 
i) and ii) for qMq and aM^ with those coming 
from a double-exponential fit of the correlator ((2.5(1 
in the time interval from t* to tf. 

Note that one could in principle consider the modified 
propagator G{t) divided by G^{t), because such a ratio 
may directly provide the mass difference a{M^ — Mq). 
However the leading exponential decays in G{t) and {t) 
occur in different time intervals, as we will illustrate later, 
causing a mismatch between the plateau of the effective 
mass of the numerator and the one of the denominator. 

In the next Section we apply our method to the extrac- 
tion of the first excited state from the correlator 1(2.5(1 . 
and we show that our extracted state can be reasonably 
identified with the Roper resonance^. 

3. LATTICE STUDY OF THE ROPER MASS 

We have adopted Wilson fermions with non- 
perturbative Clover improvement and we have gen- 
erated 450 quenched gauge configurations at /3 = 
6.20 in a lattice volume of 24^ x 64. The hop- 
ping parameter is assigned the following values: k G 
{0.1352,0.1347,0.1342,0.1337,0.1332}. The analysis of 
pseudoscalar (PS) and vector meson masses yields a crit- 
ical hopping parameter equal to kc — 0.13588(2) and an 



^ We mention that our method has been successfully applied 
also to the analysis of the two-point correlator of pseudo-scalar 
mesons. It turns out that the extracted first excited state can be 
identified with 7r(1300) resonance. However, because of both pa- 
per's length limitations and the larger phenomenological interest 
of the Af'(1440) resonance, we will present in this letter only the 
baryonic results. 



inverse lattice spacing given by = 2.55(3) GeV. Note 
that our lattice spacing is significantly smaller than the 
ones adopted in Refs. 0,0, aaMLMuJ- The statistical 
errors for all the quantities considered in this work are 
evaluated using the jacknife grouping procedure 18] . The 
values of the bare quark mass niq = (1/fc — l/kc)/2a and 
of the PS meson mass are: nig e {47,82,117,153,189} 
MeV and Mps G {0.513, 0.686, 0.831, 0.960, 1.08} GeV. 



A. Determination of the positive-parity masses 

The masses of the ground and first excited positive- 
parity states can be obtained from the correlators G{t) 
[Eq. (jTsJ] and G{t) [Eq. The ground state mass is 

extracted in a standard way by determining the plateau 
of the effective mass function corresponding to the cor- 
relator G{t) 

aM, f f it)= In {G{t)/ Git + I)} . (3.1) 

Such a mass is identified with that of the nucleon, and in 
the following indicated with aMjy. 

The sum of the masses aMpf + aAin' of the ground 
and first excited states is correspondingly obtained from 
the plateau of the effective mass function related to the 
correlator G{t) 

aMeff{t)=\n\d{t)/G{t + l)^ . (3.2) 

The quality of the plateaux for aM^f f [t) and aMe/ / {t) 
is illustrated in Fig. 

In Tabled we report the mass values aM^r and aM]^> 
obtained from the plateaux of Eqs. 1(3.1(1 and 13.2(1 . re- 
spectively, at the different values of the hopping param- 
eter used in the simulations. The time regions of the 
plateaux are reported in Table ^ and explicitly shown 
in Fig. n by the solid lines. Notice that the plateaux 
of Eq. 1(3. 2() are of comparable extension with respect to 
those of the ground state effective mass, Eq ((3.1(1 . 



k 


ti - t*f 
ti - tf 


aMjv' 
aMjv 


Mjv' /Mjv 


0.1332 


10 - 17 
20 - 30 


1.189 ±0.0181 
0.7534 ± 0.0027 


1.578 ±0.025 


0.1337 


10 - 17 
20 - 30 


1.132 ±0.022 
0.6911 ± 0.0033 


1.638 ± 0.034 


0.1342 


10 - 17 
20 - 30 


1.067 ±0.031 
0.6262 ± 0.0042 


1.705 ±0.052 


0.1347 


10 - 17 
20 - 28 


0.969 ± 0.039 
0.5561 ± 0.0058 


1.743 ± 0.075 


0.1352 


10 - 17 
20 - 28 


0.825 ± 0.053 
0.4832 ± 0.0054 


1.71 ±0.11 



Table 1: Effective masses oMn and aMj^i at the different val- 
ues of k used in the simulations. The second column reports 
the time extensions of the plateaux of Eqs. (13.111 and 13. 2t . 
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(a) k= 0.1 332 



5 10 15 20 25 30 
t 



(c) k = 0.1 342 



5 10 15 20 25 30 
t 




Fig. 1: Plateaux of the effective mass func- 
tions 13. H (diamonds) and 13.21 (dots) for k £ 
{0.1332,0.1337,0.1342,0.1347}. The solid lines are the 
time extensions considered for the plateaux of Eqs. 13.1|l and 



To check the consistency of our results, we have then 
performed a double-exponential fit of the correlator (|2.5|) 
in the time coordinate, limiting the fitting region to the 
one including the plateaux of the two effective mass anal- 
yses, i.e. [t*,tf]. In this fit we have considered two op- 
tions for the nucleon mass aMpf-. a) a free parameter to 
be fixed by the fitting procedure; b) a parameter fixed 
at the value obtained from the effective mass analysis of 
Eq. H3.1|l [see the third column of Table The two op- 
tions provide consistent results for both aMff and aAI^' , 
but the uncertainties are definitely smaller in the case b). 
Our results for aM^i, corresponding to the latter op- 
tion, are reported in Table [21 The consistency between 
the results of the effective mass analysis and the double- 
exponential fit can be appreciated from Fig. [21 where it 
can be seen that the more precise results are those ob- 
tained from the effective mass analysis. 



h 


ti — tf 




0.1332 


10 — 30 


1.154 ± 0.033 


0.1337 


10 - 30 


1.111 ±0.047 


0.1342 


10 - 30 


1.076 ± 0.074 


0.1347 


10 - 28 


1.02 ±0.11 


0.1352 


10 - 28 





Table 2: Results for the mass of the first excited positive- 
parity state obtained from a double-exponential fit of the 
correlator Ij2.5|l in the time region [tl , t /] which includes the 
plateaux of the effective masses (lO and In this fit 

the nucleon mass aMjv is constrained to be the one obtained 
from the effective mass analysis 13.111 . 




0.00 



0.05 



0.10 
(aM f 



0.15 



0.20 



Fig. 2; Comparison of the mass of the first excited positive- 
parity state, aMf^i, obtained for different quark masses us- 
ing the effective mass method (full dots) and the double- 
exponential fit (open squares). 



In Table [3 we report the values obtained for aMjy* 
from the effective mass analysis of the correlator 
together with the corresponding time intervals for the 
plateaux. The latter are shown in Fig. [31 for a couple of 
values of the quark mass. It can be seen that the ground- 
state signal is of good quality, but with respect to the case 
of positive-parity the plateaux are narrower and shifted 
to lower values of the lattice time, mainly because of the 
higher mass of the negative-parity ground state. 

In Fig. ^ our results for aM^* are plotted versus the 



B. Determination of the negative-parity mass Mjv* 

The analysis of the negative-parity ground state of the 
correlator (|2.5|) strictly follows that of the positive-parity 
one, described in the previous subsection. One simply 
has to consider t sufficiently large so that Eq. (|2.5|l can 
be replaced by 

G{t) = -J2zr e-"''-^^-'^ (3.3) 

i>0 



k 


(T-t), - {T-t)f 


aMN' 


0.1332 


14 - 19 


0.926 ± 0.009 


0.1337 


14 - 19 


0.860 ±0.011 


0.1342 


14 - 19 


0.791 ±0.015 


0.1347 


14 - 19 


0.717 ±0.023 


0.1352 


12 - 17 


0.697 ± 0.026 



Table 3: Results for the negative-parity ground-state mass 
aMiv* obtained through the effective mass analysis 13. H . The 
time intervals of the plateaux are also reported. 
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(a) k = 0,1 332 




5 10 15 20 5 10 15 20 

T-t T-t 

Fig. 3: Plateaux of the efFective mass function l|3.1^ for the 
(negative-parity) correlator 13.31 1 for k — 0.1332 (a) and k = 
0.1347 (b). 



squared PS mass [aMps)'^ ■ As a consistency check, we 
have performed a single-exponential fit of the correlator 
(|3.3|) in the time interval corresponding to the plateaux 
reported in Table |2| The results are shown in Fig. ^ as 
open squares. 
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1; 
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m 



0.00 0.05 0.10 0.15 0.20 

(aM f 

Fig. 4: Results obtained for uMn* from the effective mass 
analysis (full dots) and from a single-exponential fit in the 
time interval (open squares) versus the square of the PS mass. 
The sum of the masses uMm + aMps (open triangles) is also 
plotted for comparison. 

Our lattice simulations can be directly compared with 
those of Ref. [19], since the only difference is the use of 
smeared interpolating operators in the latter. The two 
simulations nicely agree in the results for both aM^v and 

In Fig.^the mass aA/jv + aMps of the (S-wave) scat- 
tering state N + PS at zero hadron momenta, is reported 
for comparison. Note that such a scattering state stays 
higher in energy with respect to the resonance state for 
the quark masses used in the present simulations, except 
for the lightest mass where they nearly coincide (cf. also 
Ref. 01). In principle the mass of the scattering state 
can be extracted using the modified correlator (|2.7|l con- 
structed as in Eq. H2.6|) from the negative-parity propa- 
gator H3.3|l . However in practice, due to the occurrence 



of the plateaux of the negative-parity ground-state at 
quite earlier times compared with the case of the positive- 
parity one (see Tables ^ and , as well as to a smaller 
mass gap and to larger statistical fluctuations, we find no 
clear signal of the scattering state N+PS in the modified 
correlator. 

Though it is beyond the aim of the present study, we 
want to briefly comment on the systematic errors of our 
simulations before closing the subsection. 

Our calculations have been carried out at a fixed value 
of the lattice spacing a. The use of the non-perturbative 
Clover improvement guarantees that discretization errors 
of order 0{a) are absent in our extracted masses. Nev- 
ertheless the closeness to the continuum limit strongly 
depends on the condition aM <ti 1, which is better ful- 
filled by our results for the positive-parity ground-state 
(the nucleon). In case of aM^* and uMn' the discretiza- 
tion error should be investigated by performing simula- 
tions at different values of the lattice spacing. In this 
respect we point out again that the value of a employed 
in our study is much lower than those used in most of 
the works appeared in the literature on the subject, like 
Refs. II li Id 111 Mill. 

As shown in finite volume effects can affect the 
extracted mass of excited states. The estimates made in 
indicate that at the volume and the quark masses used 
in our lattice simulations one can expect a w 10 20% 
reduction of the resonance masses. 

Finally, the error due to the quenched approximation 
can be estimated only by removing such approximation. 
Available results in the literature suggest that the 
quenching error on the ground-sate baryon spectrum is 
not large and expected within w 10%. However, very 
little is known about the impact of the quenched approx- 
imation on the excited baryon spectrum. 



C. Extrapolation to the physical point 

All our lattice results for aM^, qMn* and aM^i are 
shown in Fig. O versus the square of the PS mass. Up to 
the lightest quark mass used in our simulations there is no 
level switching between the negative-parity ground-state 
and the positive-parity first excited state, the latter lying 
above the former at variance with the physical spectrum. 

The extrapolation to the physical point should be 
guided by chiral perturbation theory. However, the 
quark masses used in our simulations are quite far from 
the chiral regime and indeed we expect that the results 
shown in Fig. [S] do not contain sizable effects from chi- 
ral logs. In order to identify the masses found on the 
lattice with the resonances of the physical spectrum, 
one should perform simulations at lower quark masses. 
For a naive correspondence we will limit ourselves to 
adopt an analytical quark mass dependence of the form 
aM = A + B{aMpsf + C{aMpsY, which is repre- 
sented by the lines in Fig. [S] for the various states. The 
positive- and negative-parity ground states can be rea- 
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^ PS 



0.15 



0.20 



Fig. 5: Plot of aA'/jv (diamonds), aM^* (squares) and aMis[i 
(dots) versus the square of the PS mass. The dotted, dashed 
and sohd hues are analytical fits of the form A + B(aMpsY + 
C{aMps)'^ ■ The markers at the physical point correspond in 
lattice units to the physical masses of the nucleon N(938), the 
5'ii(1535) and iV'(1440) resonances. 



sonably identified with nucleon A^(938) and the S'n (1535) 
resonance. As for the positive-parity first excited state, 
the quark mass dependence of our extracted masses is 
consistent with the identification with the Roper reso- 
nance Af'(1440), though a possible interpretation as the 
higher nucleon resonance Ai"(1710) cannot be completely 
ruled out. Such a possible extrapolation ambiguity is 
not addressed in the literature 8, 9, 13, 14, 15], and a 
clear-cut solution again relies in performing lattice sim- 
ulations at lower quark masses. The latter would also 
help clarifying at which quark masses the level switching 



between the orbital and radial excitations of the nucleon 
occurs. Basing on our crude extrapolations such a cross- 
ing is tentatively expected to appear close to the physical 
point. 



4. CONCLUSIONS 

We have presented a new method for extracting ex- 
cited states from a single two-point correlation function 
calculated on the lattice. We have defined a new corre- 
lator which is simply an appropriate combination of the 
correlation function evaluated at different time slices. In 
this way the leading exponential decay (ground state) 
is exactly subtracted and one can access the first ex- 
cited state. The mass of this state can be extracted from 
the plateau of the effective mass of the new correlator. 
Therefore our method allows to establish unambiguously 
the time interval where the ground and the first excited 
states dominate the two-point correlation function. 

Our method has been applied to a quenched lattice 
study (volume = 24^ x 64, /3 = 6.2, a-^ = 2.55 GeV) 
of the first excited state of the nucleon using the inter- 
polating operator O — eabc[u^Cl5db]uc- The results are 
consistent with the identification of our extracted excited 
state with the Roper resonance iV'(1440). The switching 
of the level ordering with respect to the negative-parity 
partner of the nucleon, iV*(1535), is not seen at the simu- 
lated quark masses and, basing on crude extrapolations, 
is tentatively expected to occur close to the physical 
point. 
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